Influence of statistical sequential decay on isoscaling and symmetry energy coefficient 

in the GEMINI simulation 
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Extensive calculations on isoscaling behavior with the sequential-decay model GEMINI are per- 
formed for the mediate-heavy nuclei in the mass range A = 60-120 at excitation energies up to 3 
MeV per nucleon. The comparison between the products after the first-step decay and the ones 
after entire-steps decay demonstrates that there exists strong sequential decay effect on the final 
isoscaling parameters and the apparent temperature. Results show that the apparent symmetry 
energy coefficient jap P does not reflect the initial symmetry energy coefficient C sym embedded in 
the mass calculation in the present GEMINI model. 
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One of main goals of the isospin physics is to determine 
the isospin dependence of the in-medium nuclear effective 
interactions and the equation of state (EOS) of isospin 
asymmetric nuclear matter or finite nuclei, particularly 
its isospin-dcpcndent term, i.e., the density dependence 
of the nuclear symmetry energy. Knowledge of nuclear 
symmetry energy is essential for understanding not only 
many problems in nuclear physics, such as the dynam- 
ics of heavy-ion collisions induced by radioactive beams 
and the structure of exotic nuclei, but also a number of 
important issues in astrophysics, such as the supernova 
simulation and neutron star models, which require inputs 
of the nuclear equation of state at extreme values of den- 
sity and asymmetry [1, 2]. Recently impressive progress 
has been made both experimentally and theoretically, a 
number of earlier reviews on isospin physics with heavy- 
ion reactions can be found in several references [3-5] . 

Symmetry energy could be extracted from heavy-ion 
collision using the isoscaling approach [6-15]. Isoscaling 
law means that the ratio of isotope yields R2i(N, Z) = 
Y2(N, Z)/Y\(N, Z), from two similar reactions, denoted 
as reaction 1 and 2, which are different only in their 
isospin asymmetry, is found to exhibit an exponential 
relationship as a function of the neutron number N and 
proton number Z [6], i. e. 



R 21 {N, Z) 



Y 2 (N, Z) 
Yi(N, Z) 



Cexp(a7V + /3Z), (1) 



where Y 2 (N, Z) and Y\ (N, Z) are the fragment yields 
from the neutron-rich and the neutron-deficient reaction, 
respectively, C is an overall normalization factor, and a 
and {3 are fitted parameters. The isoscaling parameter 
a is related to the symmetry energy coefficient C sym of 
EOS in microcanonical and canonical frames by following 
relation [6], 
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where Z s \, Z s2l N„i, N s2 , A s \, A s2 are the charge num- 
ber, neutron number and mass numbers of the sources 
from the two systems, T is their temperature and C sym 
is the symmetry energy coefficient. This relation has 
also been evidenced in other model frameworks [7-15]. A 
great deal of effort has been devoted to investigate the nu- 
clear symmetry energy and its density/temperature de- 
pendences [3, 5, 16-19]. 

Ideally, primary fragments should be detected right af- 
ter emission in order to extract information about the 
collisions, and Eq. (2) is derived based on the primary 
reaction products bypassing secondary decays. However, 
the detected experimental data are for cold products af- 
ter the secondary decays from hot products. Isoscaling 
has also been reasonably reproduced in the sequential 
decay codes [15, 20]. However there are still arguments 
on the sequential decay effect on isoscaling, some models 
show that the effect from sequential decays on isoscaling 
is negligible, but some efforts show that sequential decay 
affects on the isoscaling parameters, and then distort the 
extraction of symmetry energy coefficient C sym [21, 22]. 
There are some issues still keep unsolved or unclear, such 
as the sequential decay effect on isoscaling parameters, 
derived apparent temperature T r from the experimental 
measurement and the isospin evolution of the decaying 
sources, these factors affect the extraction of symmetry 
energy coefficient C sym . 

The statistical GEMINI model [23] calculates the de- 
cay of compound nuclei by modes of sequential binary 
decays. The model employs a Monte Carlo technique 
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to follow the decay chains of individual compound nu- 
clei through sequential binary decays until the resulting 
products arc unable to undergo further decay. GEM- 
INI has been widely used to simulate the hot equilibrium 
source de-excitation, or as an "afterburner" code to ana- 
lyze the hot fragments decay after dynamical simulation 
[24-26]. Isoscaling has been investigated by statistical 
sequential secondary decay code GEMINI [23], in which 
only the first-step sequential decay was simulated and 
Eq. (2) was confirmed for the fragments which are de- 
cayed directly from the initial sources [15]. In the present 
work, we investigate the entirely decayed fragments from 
excited sources, comparing with the only first-decay frag- 
ments from the same source. The influences of sequential 
decays on isoscaling parameters a, ft, and the apparent 
temperature T r are discussed, and the apparent symme- 
try energy coefficient 7 app is extracted. 

The detailed description of GEMINI code can be found 
in [15, 23], the same configuration and parameters of the 
GEMINI code were adopted as in Ref. [15]. Several 
pairs of equilibrated sources are considered at various 
initial excitation energy E ex ~1.0, 1.4, 2.0, 2.4, and 3.0 
MeV/nuclcon. We selected source pairs with the same 
proton number Z s but different mass number A s to sys- 
tematically study the isoscaling behavior. In this case, 
possible effects of different magnitudes of Coulomb inter- 
action on isotopic distributions are avoided. The equi- 
librated source pairs are chosen in different mass region 
and system isospin asymmetry N/Z. Two groups of the 
source pairs have been used: (1)Z S =50 with ^4 S =100, 
105, 110, and 115, respectively; (2)Z S =30 with A,=60, 
63, 66, and 69, respectively Following literature the 
index "2" denotes more neutron-rich system as widely 
used in convention, and index "1" denotes the the more 
neutron-deficient system. In our previous work [15], the 
statistical decay stops after one particle emitted from the 
source, that was called the "first-step" decay in this pa- 
per, which is a simple picture and the decay procedure 
can be expressed definitely and clearly, isoscaling has 
been confirmed for the first-step decay products in de- 
tail, and the reasonability of extracting symmetry energy 
coefficient C sym from the simulation results via experi- 
mental analysis technique. But the first-step decay was 
not the real case, experiments measure the final prod- 
ucts after multi-step decays until no fragments produced 
or gamma rays emitted, which was called "entire-steps" 
decay in this paper. 

Isoscaling is analyzed from the emitted light fragments 
in above both cases, namely first-step decay only and 
entire-steps decay chains for all the simulated systems. 
As an example shown in Fig. 1 the comparison of the 
isoscaling parameters a and ft is plotted as a function 
of emitted light fragment proton number Z and neutron 
number N, between the first-step decay and entire-steps 
decay fragments for the source pairs of Z s = 50. As we 
can see, isoscaling parameters a and ft are essentially in- 
dependent on Z or N of fragments, for both first-step 
and entire-steps decay fragments, especially for Z < 4 or 
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FIG. 1: (Color online) Comparisons of isoscaling parameters a 
(positive values) or /3 (negative values) as a function of Z or N 
from source pairs of Z s = 50 at excitation energies E ex = 2.4 
MeV/nuclcon. All solid symbols represent the results for only 
the first-step secondary decay products, open symbols for entire- 
step secondary decay products. Symbols in the figure correspond 
to Y A s =115 /Ya s =ioo (up-triangles), Y As=no /Y As= i 00 (circles), 
Y As= no /Y As= ios (down-triangles), respectively. 

N < 5. It has been evidenced that in the first-step decay 
case, the probability of producing a cluster with a given 
Z and A at T depends exponentially on the free energy 
of that cluster, F(Z,A,T), in GEMINI simulation, the 
cluster free energies depend on the strength of the sym- 
metry term in the liquid-drop energy through the Eq. 
(2) with C sym ~ 24 MeV[15]. If the entire-steps sequen- 
tial secondary decay is included, the source isospin N/Z 
and temperature T as well as isoscaling parameters vary 
after each step of decay. Similarly, in previous study, 
time dependence of isoscaling parameters have been dis- 
cussed in molecular dynamics model [14], indicating that 
the final values of these parameters could be related to 
the last part of the reaction where the fragments fin- 
ish cooling by particle evaporation. In present GEMINI 
model, fragment yields are strongly expected degraded 
after entire-steps decay, so do the isotopic yield ratios. It 
is not surprising that in Fig. 1 that parameters a and ft 
extracted from isotopic yield ratios of the final emitted 
light fragments show discrepancy, especially for the inter- 
mediate mass fragments like Z > 5 in some cases, since 
those heavier fragments experience strong multi-step de- 
cay and feeding-down effects. Finally, a and ft change a 
lot comparing with the only first-step decay case, and the 
fluctuation increases in the entire-steps decay case. Av- 
erage a and ft values over fragments Z and N are used 
in following discussions to observe the overall property. 

In Fig. 2 we show the comparison of isoscaling pa- 
rameters a and ft as a function of excitation energy of 
sources from different source pairs. If the entire-steps 
decay chains are included, which are depicted by the 
open symbols in Fig. 2, a and \ft\ values show signifi- 
cant decrease. This reduction is about 20% in average, 
consistent with the result in Rcf. [20, 21] for the C sym f» 
25MeV case, but the excitation energy dependent trend 
of a and ft does not change. 
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FIG. 2: (Color online) Comparisons of isoscaling parameters a 
(left panel) and (right panel) as a function of the source excitation 
energy from source pairs of Z s =30. All solid symbols represent only 
the first-step secondary decay products, open symbols are entire 
steps secondary decay products. Symbols in the figure correspond 
to Ya s =69 /Ya 3 =60 (down-triangles), Ya s =69 /Yk s =63 (diamonds), 
Y A e =m/ Y A s =&& (left-triangles). 
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As already discussed, the isoscaling parameter a is re- 
lated to the symmetry energy coefficient C sym of the 
nuclear binding energy through Eq. (2) and (3), this 
relation provides a direct link between the measurable 
quantities and the nuclear symmetry energy coefficient. 
It should be noticed that the parameters a and (3 refer 
to the hot primary fragments, which have to undergo se- 
quential decays into the cold fragments. It was assumed 
that the secondary decay on the yield of a specific isotope 
is similar for the two reactions, thus the effect of the se- 
quential decays on R 2 i(N, Z) is small and that R,2i{N, Z) 
reflects the properties of the primary source. In our 
present study, the first-step statistical sequential decay 
process stems from a fixed initial source with definite 
excitation energy (temperature) and isospin asymmetry, 
and it has been verified by Eq. (2) and Eq. (3), to reflect 
the link between the measurable quantity R,2\{N, Z) and 
the symmetry energy coefficient C sym . 

It has been found in Fig. f , the isoscaling behavior still 
presents after considering the entire-steps decay chains. 
But in this case, the isoscaling parameters a and f3 de- 
crease comparing with the only first-step decay case as 
seen in Fig. 2. In the statistical sequential decay, the 
source temperature T and isospin asymmetry N/Z also 
change after each step of the decays, thus the parameters 
T and A{Z/A) 2 S (or A(N/A) 2 ) varies during the sequen- 
tial decay process, where many intermediate sources are 
different from the initial source. 

To explore the validity of Eq. (2) and (3) in the entire- 
steps statistical sequential decay, we plot aT as a func- 
tion of A(Z/A)l and /3T as a function of A(N/A) 2 S in 
Fig. 3. For the first-step only decay, data points de- 
picted by the solid symbols, with using the initial source 
temperature Tj which was calculated by input excitation 
energy as shown in Table I and isoscaling parameters af- 
ter the first step decay. The linear fit (the solid line) gives 



FIG. 3: (Color online) aT (/3T) as a function of the source isospin 
difference A(Z/A) 2 S (A(N/A) 2 S ) from source pairs of Z s = 30. All 
solid symbols represent only the first-step secondary decay case, 
open symbols without /with cross represent the entire-steps sec- 
ondary decay products with Ti and T r , respectively. Symbols in 
the figure correspond to excitation energies E ex /A = 1.0 (squares), 
1.4 (circles), 2.0 (up- triangles), 2.4 (down-triangles) and 3.0 MeV 
(diamonds), respectively. The solid, dash and dot lines are the 
linear fitting for the above three cases which gives the apparent 
symmetry energy coefficient fapp = 24.2, 19.65 and 15.84 MeV, 
respectively. 

a symmetry energy coefficient C sym = 24.2 ± 0.3MeV. 
To investigate the decay effect on symmetry energy coef- 
ficient, data points from the entire-steps decays are also 
plotted in Fig. 3 as shown by the open symbols. As we 
have mentioned for the entire-steps decay chains, there 
are many intermediate-state sources with different tem- 
perature T and isospin asymmetry Z/A or N/A. Nev- 
ertheless, the final isoscaling parameters a and /3 still 
show similar rules as the first-step decay only case, i.e. 
aT (0T) and A(Z/A) 2 S (A(Z/A) 2 S ) can be still fitted by 
another linear function, namely 




.(5) 

which gives the apparent symmetry energy coefficient 
lapp = 19.65 ± 0.25 MeV if the same T^ are used. 

In fact, for the source decays with entire steps, the 
intermediate-state sources have different isospin asymme- 
try N/Z ranging from initial isospin asymmetry to the 
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TABLE I: Initial temperature (Ti) and final-state apparent 
temperature (T r ) in different systems (Z 3 —30, A s — 60 or 63, 
66 or 69. 





Ti 


T r 


T 


T r (MeV) 


A s 


60/63 


60/63 


66/69 


66/69 


1.0 


2.8/2.9 


2.3/2.4 


2.9/2.9 


2.4/2.4 


1.4 


3.4/3.5 


2.8/2.9 


3.5/3.5 


2.9/3.0 


2.0 


4.1/4.2 


3.5/3.6 


4.2/4.2 


3.6/3.6 


2.4 


4.5/4.6 


3.8/3.9 


4.6/4.7 


3.9/4.0 


3.0 


5.1/5.1 


4.3/4.4 


5.2/5.3 


4.4/4.5 



stable line or the evaporation attract line [27] and dif- 
ferent temperature T ranging from initial temperature 
to zero. In this case, principally both T and A(Z/A) 2 S 
(A(Z/A)1) need to be corrected to reflect the inter- 
mediate sources. In the simulation, the initial source 
temperature T can be calculated, and the intermediate 
source tracing the sequential decay chains can also be 
performed. From an experimental point of view, temper- 
ature and isospin asymmetry of the intermediate source 
can be extracted from evaporation products which re- 
flects the whole decay chains. 

Traditionally, temperature can be extracted from the 
measurements of spectral slopes or double isotopic ratios 
at lower energies [28, 29]. In the present work, initial 
temperature Ti are calculated directly in the GEMINI 
code by the input excitation energy [15], and the final- 
state temperature T r can be obtained by the neutron and 
proton spectra fitting when the entire-steps decay chains 
are included in the GEMINI calculation. The results are 
displayed in TABLE. I. 

When the temperature T r is used in Eq. (4) and (5) 
to fit the linear slope parameter 7 app , it leads to the re- 



duction of the parameter 7 app as shown in Fig. 3 (short 
dashed line). Its slope gives an apparent symmetry en- 
ergy coefficient -f app = 15.84 ± 0.18 MeV, which is one- 
third reduction comparing with the only first-step decay 
case. In this context, we should be careful to use the ap- 
parent symmetry energy derived directly from the final 
fragments which could be distorted due to the multi-step 
sequential decays. Of course, the present results are spe- 
cific for the use of GEMINI to describe the secondary 
decay, i.e. they may depend on the details of the sequen- 
tial decay code. 

In summary, we performed the isoscaling analysis for 
both light fragments from only the first-step decay and 
the entire-steps decay chains with GEMINI code, it is 
found that isoscaling can still be observed and the Eq. 
(4) and (5) which are used to extract the symmetry en- 
ergy coefficient also work after the entire-steps decay is 
taken into account. However, the statistical sequential 
decay leads to the decreasing of isoscaling parameters 
a and ft as well as temperature. Therefore, the reduced 
(apparent) source temperature together with the reduced 
isoscaling parameters leads to a smaller symmetry en- 
ergy parameter 7 app in comparison with the initial sym- 
metry energy coefficient C sym which is constrained from 
the only first-step statistical decay calculation. From the 
present GEMINI model calculations, we shall carefully 
consider the multi-step sequential decay effect on the ex- 
traction of the symmetry energy coefficient via the final 
cold products. 
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